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Abstract 

In this paper we develop a dressing method for constructing and solving some classes 
of matrix quasi-linear Partial Differential Equations (PDEs) in arbitrary dimensions. This 
method is based on a homogeneous integral equation with a nontrivial kernel, which allows 
one to reduce the nonlinear PDEs to systems of non-differential (algebraic or transcen- 
dental) equations for the unknown fields. In the simplest examples, the above dressing 
scheme captures matrix equations integrated by the characteristics method and nonlinear 
PDEs associated with matrix Hopf-Cole transformations. 

1 Introduction 

Since the pioneering work [1] on the Korteweg - de Vries equation [2], completely integrable 
nonlinear Partial Differential Equations (PDEs) have been intensively studied during the last 
decades, and large classes of integrable PDEs have been found, like the so called S-integrable 
systems (or soliton equations, whose integration scheme involves the solution of a linear integral 
equation) [3, 4], and the so-called C-integrable equations (integrated by simpler transforma- 
tions, like the Hopf-Cole transformation for the Burgers equation) [5]. Much effort has been 
devoted to the study of direct techniques to construct and solve nonlinear PDEs. One of the 
most powerful of such techniques is the dressing method, originally developed for (1+1) and 
(2+l)-dimensional S'-integrable models [6, 7, 8, 9] (see also [10]). Multidimensional general- 
izations of it have also been developed [11, 12, 13, 14], allowing to integrate special classes of 
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higher dimensional nonlinear PDEs. Nonlinear dressing methods for PDEs in arbitrary dimen- 
sions associated with vector fields are also known [15, 16]. All the above dressing formalisms 
are based on the hypothesis that, for given spectral data, the spectral function can be uniquely 
constructed from the relevant integral equation; i.e., the kernel of the corresponding integral 
operator is empty. 

Recently, a new version of the dressing method has appeared, based on integral operators 
with nontrivial kernel [17]. This assumption removes most of the restrictions on the dimen- 
sionality of the space of analytic solutions of the constructed PDEs, and the solutions of the 
n-dimensional PDEs constructed in [17] are parametrized by arbitrary spectral functions of 
(n — 2) variables (full integrability would be achieved if the variables were (n— 1)). The struc- 
ture of the nonlinear PDEs has been simplified and the space of analytic solutions has been 
enriched in [18], where multidimensional dressing operators [19], which are not effective in the 
classical case, have been used. The solutions of the n-dimensional nonlinear PDEs constructed 
there are parametrized by arbitrary spectral functions of (n— 1) variables, but full integrability 
has not been achieved even there, since these spectral functions are not in the right number. 

In this paper we develop another variant of the dressing method, based on a homogeneous 
integral equation with nontrivial kernel, allowing one to reduce certain classes of nonlinear 
PDEs in arbitrary dimensions to systems of non-differential (algebraic and/or transcendental) 
equations, similarly to the method of characteristics [20]. The nonlinear PDEs in arbitrary 
dimensions isolated by this method are built in terms of "dressed" first order operators of the 
type (4) below. To increase the dimensionality of these PDEs and of their space of analytic 
solutions, we use multidimensional differential operators first introduced in [18]. 

Below is the list of the nonlinear PDEs which will be derived and solved in this paper. All 
equations are Q xQ matrix equations (Q G N + ), unless differently specified. Hereafter we write 
superscripts inside of parenthesis in order to distinguish them from the power notation. 

1. The first order matrix equation 

N 

w tl +J2 w *iP U) ( w ) = [w,Tp {0 \w)} (1) 

3=1 

will be derived in Sec. 3. Here T is any constant matrix, p^(w) are arbitrary matrix functions 
representable as positive power series of w, and N is any integer. This matrix equation can 
also be integrated [21] by the method of characteristics. 

2. A minor modification of the dressing method for eqs.(l) allows one to construct a class of 
second order matrix systems (see Sec. 5) which degenerate, for the simplest choice of the arbi- 
trary functions, to the matrix Burgers equation, linearizable by the Hopf-Cole transformation. 
A typical example is given by the second order matrix system: 

C 2 (w) = [w,S(v 2 - d(v))v], , , 

C 2 (v) + SC\{v) = SC^v 2 ) + [SCi(v),v] + [v, Sv% {l} 

subjected to the constraint 

Ci(w) = [w,v], (3) 

where S is any constant diagonal matrix, the differential operators C m , m = 1,2, are defined 
as follows: 

N 

cum) = f tm (x) +^/*>)p (mj) w, (4) 
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and p^ m ^{w) are arbitrary matrix functions representable as positive power series of w. 

We remark that many equations of Mathematical Physics appear as first order quasilinear 
PDEs in multidimensions. Therefore it is important to develop efficient methods to isolate and 
solve, among these equations, integrable and/or partially integrable cases. 

The matrix equation (1), a remarkable example of integrable system of PDEs in arbitrary 
dimensions, is the natural matrix generalization of the classical examples of physically relevant 
scalar equations integrable by the method of characteristics. A systematic study of the matrix 
reductions of (1), with the goal of isolating physically relevant cases, will be the subject of a 
subsequent paper. Here we only remark that, in the 2x2 matrix reduction 

w=(* (5) 

equation (1), with N = 1 and p^ = 0, coincides [21] with the gas-dynamics equations [20] 

V\ t + V\V\ Xl + V2V2 Xl = 0, , g * 
V2t + Wlxi + V1"2 X1 = 0. 

We also remark that other vector generalizations of scalar equations integrable by the method 
of characteristics have been studied during the last twenty years, by different methods, in a set 
of papers [22, 23, 24]. 

The system of second order equations (2), constructed in terms of the two first order multidi- 
mensional operators (4), presents mathematical features in common with equations integrable 
by the method of characteristics, and with equations integrable by the Hopf-Cole transfor- 
mation; therefore it is conceptually similar to the equations, introduced in [14], presenting 
mathematical features of both S - integrable and C - integrable PDEs. Equations like (2), 
for which integrability properties of different type merge together, are interesting prototype 
examples in the theory of integrable systems. 

The paper is organized as follows. In Sec. 2, after a brief review of the classical dressing 
method (in Sec. 2.1) and of the novel dressing features contained in [17] (in Sect. 2. 2), we describe 
the main features of the general dressing algorithm of this paper. Eqs.(l) will be derived in 
Sec. 3. In Sec. 4 we use the dressing algorithm to reduce Eqs.(l) to a system of non-differential 
equations characterizing the general solution of (1); in particular, we discuss the Cauchy prob- 
lem. In Sec. 5 we derive eqs.(2) and the associated general class, discussing their solution space. 
In Sec. 6 we make some remarks on the overdetermined system of linear PDEs for the associated 
spectral function. Conclusions are presented in Sec. 7. 

2 Basic novelties of the dressing method 

To emphasise all the significant novelties and features of our algorithm, in comparison with both 
the classical dressing method and the new dressing method developed in [17], first, we give a 
brief review of one of the versions of the classical dressing method for the (2+l)-dimensional N- 
wave equation (Sec.2.1), and, second, we describe the modifications introduced in [17] (Sec. 2. 2). 
After such a preliminary overview, we explain the novelties of the algorithm introduced in this 
paper, (Sec. 2. 3). 
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2.1 Brief review of the classical version of the dressing method, 
and the (2+l)-dimensional iV-wave equation 

The starting point of one of the versions of the classical dressing method is the linear integral 
equation 

$(A; x) = J *(A, /x; x)U(n\ x)dtt(fi) = VU(\; x), (7) 

in the spectral variables A = (Ai, . . . , AdimA), A* = (a*i 5 • • • , fi><Km\), f° r the unknown matrix 
function U. The given matrix functions $ and ^ are defined by some extra conditions, which 
fix their dependence on an additional vector parameter x — (x±, . . . , Xdimx), whose components 
are the independent variables of the associated nonlinear PDEs. Q is some largely arbitrary 
scalar measure in the A^-space. Apart from Q, all the functions appearing in this paper are 
Q x Q matrix functions. 

We remark that no a priory assumption is made in (7) on the dependence of \I/ on A (this 
general starting point has been used, for instance, in [25] and in [14]), to keep the structure 
of ^ as much general as possible. Indeed, although in most of the cases such a dependence is 
described by a Cauchy kernel, an indication that equation (7) is the manifestation of Riemann- 
Hilbert and/or d analyticity problems, there are examples (see [25] and [14]) in which more 
general representations appear, indicating that the above analyticity problems could be a too 
restrictive starting points. 



2.1.1 Derivation of iV-wave system. 

The basic assumption underlying all the known classical dressing procedures available in the 
literature, is that the operator \l/ in (7) is uniquely invertible; i.e., that 

dimker^ = 0. (8) 

The x-dependence is introduced by the matrix equations 

* Xi (A, fx;x) = $(\;x)B i C(n;x), i = 1, .., dim x, (9) 

showing that the ^-derivatives of the kernel \1/ are degenerate matrix functions, another basic 
feature of all known classical dressing algorithms, where Bi, i = l,..,dim x, are constant 
diagonal matrices, at most Q of them may be independent. Due to the above degeneracy, the 
compatibility of equations (9) leads to separate equations for $ and C: 

'I', 1>j 'I'r 1>: 0. i^j, (10) 

i>,c, /u:,, o. i^j, (ii) 

and one equation is the adjoint of the other. Without loss of generality we assume B 1 = J, 
where I is the identity matrix. 

Replacing, in equation (10), $ by ^/U, as indicated in (7), and using (9), one obtains the 
following equation: 

VLijU = 0, (12) 

where 

LijU = U Xz Bj - U Xj Bi + UBivBj - UBjvB h i,j = 1, .., dim x, i ^ j (13) 
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and 



v(x) = / C(X;x)U{X;x)dQ(X) 



(14) 



Then the property (8) implies that 



LijU(\; x) = 0, i, j = 1, .., dim x, 



i 1^ 3 



(15) 



or, explicitly: 



L 21 U = U X2 -U Xl B 2 -U[v,B 2 ] 
L 31 U = U X3 -U Xl B 3 -U[v,B 3 ] 



0, 
0, 



(16) 



having chosen j — 1, i — 2, 3. 

This is nothing but the well-known linear overdetermined system corresponding to the TV- 
wave equation in the three variables x±,x 2 ,x 3 . 

Its compatibility condition yields 



It is important to remark that 

1. eq.(17) may be derived in a different way, "saturating the parameter A" in equations (16) 
by the integral operator f dQ(X)C(X;x)-, and obtaining the nonlinear system: 



Eliminating the extra field v\ from these two equations, we get the (2 + l)-dimensional 
iV-wave system (17). 

The possibility to derive integrable systems in these two alternative ways is important, 
since, while integrable PDEs in 2+1 dimensions (or less) are characterized as the compat- 
ibility condition of a linear overdetermined system of PDEs, such a basic property seems 
to be lost in multidimensions. 

2. Each linear equation (16) is two-dimensional. 
2.1.2 Solution space 

We now consider the manifold of analytic solutions of equation (17) generated by the dressing 
procedure. The solutions of eqs.(10) and (11) can be parametrized as follows: 



[v X3 ,B 2 ] - [v X2 ,B 3 ] + B 2 v Xl B 3 - B 3 v xl B 2 - [[v, B 2 ], [v, B 3 }] = 0. 



(17) 



L 21 v - [B 2 , vi] = v X2 - v Xl B 2 - v[v, B 2 ] - [B 2 , vi] = 0, 
L 31 v - [B 3 , vx] = v X3 - v Xl B 3 - v[v, B 3 ] - [B 3 , = 0, 

written in terms of the square matrix fields v(x) and vi(x), where 



(18) 




(19) 




(21) 



(20) 
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where 

dim x 

B-x=J2 B ^ ( 22 ) 
i=i 

and where the spectral parameters A, /i, k, q are scalars. It is simple to see, from the linear limit, 
that the solution space of equation (17), generated by the dressing algorithm, is full. Indeed, 
in the linear limit: \I/ ( A, fi) ~ 5(X — fi) and U ~ <E>. Take Cq(X, q) = 5(X — q); then the solution 
v of the 3 - dimensional iV-wave system (17), which in the linear limit reads 

v{x)~ J C(X;x)${X;x)dQ(X) = J e AB ^$ (A, k)e kBx dkdtt(X), (23) 

is parametrized by the arbitrary matrix function $o(A, k) of the two scalar spectral parameters 
A, k] then its solution space is 2 dimensional, and therefore it is complete. 

2.2 Novelties of the dressing methods introduced in [17] 

In [17] we assumed that the kernel of the operator ^ is one dimensional: 

dimker^ = 1, (24) 

i.e., the solution of the homogeneous equation associated with eq.(7) is nontrivial: 

= ^H H(X;x) = U h (X;x)f(x), (25) 

where U h (X;x) is some nontrivial solution of the homogeneous equation ^/H = and f(x) is 
an arbitrary matrix function of x. Then the general solution of eq.(7) reads 

U(X;x) = U p {X;x) + U h (X;x)f(x), (26) 

where U p (X;x) is some particular solution of (7). 

As a consequence of the novel assumption (24), equation (12) implies the following equations 
for U: 

£j(X; x) = L jl U(X; x) - (L 2 iU(X; x))A j (x) = 0, j = 3, . . . , dimx, (27) 
L n U = U Xj - U X1 B 3 - U[v, Bj], j = 2, . . . , dim x, (28) 

where A^(x) are some matrix functions to be defined. We have established that, if dimker^ = 1, 
then each linear equation (27) for the spectral function U (A; x) is 3 dimensional. 

The associated nonlinear equations are obtained "saturating the parameter A" in equations 
(27) by the integral operator f dfl(X)C(X; x)-. In order to express A^{x) in terms of U and 
close the system, we introduce an external dressing function G(X;x), and the associated new 
matrix fields 

w 00 (x) = J G(X;x)U(X;x)dn(X), w j0 (x) = J G xj (X;x)U(X;x)dn(X), j > 0, (29) 
w ij (x) = J G xlx3 (X;x)U(X;x)dtt(X), i,j>0, w ij (x) = w ji (x). (30) 

Since the dimensionality of G has no formal restrictions, the above w - fields increase the 
dimensionality of the nonlinear PDEs. This can be seen, for instance, from their small field lim- 
its: w 00 (x) ~ J G(X;x)<S>(X;x)dfl(X), w j0 (x) ~ / G Xj (X; x)$(A; x)dQ(X). Nonlinear equations 
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for these fields appear after applying the integral operators f dQ(\)G(\; x)-, f d£l(\)G x .(\;x)- 
to eq.(27). 

To close the system of nonlinear PDEs, one needs (a) equations defining G(X; x) and (b) an 
additional relation between all the matrix fields, which fixes the arbitrary function f(x) in the 
solution space, see eq.(26), and may be taken in quite arbitrary form 

F(v, v\ w 00 , w l °, w ij ) = 0, i,j = 1,2,.... (31) 

Let us collect the basic novelties of the algorithm. 

1. The existence of a nontrivial kernel of the basic integral equation implies that the solutions 
constructed by the dressing depend on an arbitrary function f(x) of the coordinates; this 
fact has the following important implications. 

2. The nonlinear system of PDEs constructed by the dressing scheme possesses a distin- 
guished block structure and is underdetermined. 

3. To close the system and fix its underdeterminacy (or, equivalently, to fix f(x)), one has 
to introduce an "external and largely arbitrary" relation among the fields (see (31)). 

4. The system of PDEs depends on two types of matrix fields, those obtained "saturating 
the parameter A" of the solution U(X] x) of the linear integral equation by ingredients of 
the classical dressing method, whose dimensionality is constrained, and those obtained 
saturating A by a novel dressing function G(A; x), whose dimensionality is not constrained. 
That's why the dimensionality of the solution space, (n — 2), can be arbitrarily large. 

5. While integrable PDEs in low dimensions (2+1 or less) are the compatibility of overdeter- 
mined systems of linear spectral problems, such a feature seems to be lost for our higher 
dimensional examples. 

A prototype example of the above construction is given by the following 4 dimensional 
system of two matrix equations 

# 2 (gi, qi, ?2)^ 1 (gi, 92, 9s) = #3(51, Qi, Q2)Ba 1 (qi, q 2 , 9a) = B 4 (gi, qi, q 2 )B^ 1 (q 1 , q 2 , q 3 ) (32) 

for the three square matrix fields qi(x) , q 2 (x) , q 3 (x) , supplemented by the "largely arbitrary" 
relation 

F(qi,q 2 ,q 3 ) = (33) 

among them, where the matrix blocks Bj are defined as: 

Bj(Qi, Q2, Qs) = q 2X] - q 2xi Bj - q 2 [qi, Bj] - [Bj, q 3 ], j = 2, 3, 4, (34) 

and Bj, j = 2,3,4 are constant diagonal matrices different from the identity. In the simplest 
case, the largely arbitrary relation (33) can be chosen to be an equation defining one of the 
fields, say q 3 , to be any given function 7(2;) (in general, a generalized function), interpretable 
as an "external arbitrary forcing" : 

F: q 3 (x)= 1 (x). (35) 

The partially integrable nonlinear PDEs (32-35) possess a manifold of analytic solutions of 
dimension 2. 
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2.3 Dressing method based on homogeneous Integral equation 

Although the starting integral equation used in [17, 18] is still inhomogeneous, the use in [18] of 
first order multidimensional dressing operators allows one to enrich the space of analytic solu- 
tions of the constructed nonlinear PDEs (although full integrability is not achieved even there), 
and to simplify the structure of such nonlinear PDEs (which become differential polynomials), 
if compared to the block structure of equations like (32-35). 

In addition, the introduction of such first order dressing operators makes clear that the 
inhomogeneous term <3> of the integral equation (7) is not necessary anymore, suggesting the new 
scenario, discussed in this paper, of a dressing algorithm based on the following homogeneous 
integral equation: 

= J^(X,u;x)U(u;x)dn(u) = ^f(X,u;x)*U(u;x), (36) 
supplemented by the generalized commutation relation 



A(X, v) * /I] x) = *(A, v\ x) * A(u, fi). (37) 

In the integral equation (36), U is the unknown spectral function depending on the sin- 
gle spectral parameter A; the kernel \I> of the integral operator, the dressing function, sat- 
isfies the generalized commutation relation (37) for a proper choice of the auxiliary func- 
tions A and A. Function \1/ and, consequently, U, depend on an additional set of variables 
x = (ti, t2, ■ ■ ■ , xi, X2, ■ ■ ■ ), which are the independent variables of the associated nonlinear 
PDEs. In this paper we assume that the functions A and A do not depend on x. In general, 
all functions are Q x Q matrices. 

Following [17], we assume in this paper that the integral equation (36) possesses nontrivial 
solutions, namely, that 

dimker(^*) = d > 0. (38) 
The general solution of the homogeneous equation (36) reads 

d 

tf(A;s)=£tf«(A; *)/«(*), (39) 
i=i 

where uj^ are independent particular solutions of (36) and /W are arbitrary functions of x. It 
is convenient to introduce a "rectangular" integral operator in the following way. Let D be a 
set of points: D = . . . , l M }, and let V be a disjoint set (V fl D — 0), consisting eventually 
of continuous and discrete parts. Thus, in the function *(A, x), X e V, while /j, e V U D: 

As a consequence of this assumption, A e V U D in U(X; x): 

U(X;x) = \ U ^^ (41) 
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and the integral equation (36) reduces to the form 

M 

fi; x) * u(fi; x) + ^^ip 0i (\; x)ui(x) = 0, A G V. (42) 
i=i 

If the integral operator ip(X, fx;x)* is invertible, the solution u(X;x) is uniquely expressed in 
terms of the M arbitrary functions Ui(x), which may be identified with the functions f^(x) in 
eq.(39), and d = M. If ip(X, /x; x)* is not invertible, then d > M. As in [17], we introduce an 
external dressing function G(X; x) in the next section in order to fix f^ l \x). 

We remark that the integral equation (42) can be viewed as an inhomogeneous integral 
equation with an inhomogeneous term depending on M arbitrary functions. In the rest of the 
paper we find it more convenient to work with the homogeneous form (36). 

The rectangular structure of \& implies also different "square" structures for the integral 
operators A* and *A in (37): 

^,,) = {f >■ ^.eA (43, 



{a(A, /i), A, ji E V 
a (A), X eV, fi = l 
a n0 (/i), X = l n , fieV K J 

Qjnmi X Z n , /i l m 

where n,m = 1, . . . , M. 

To end this section, let us carry out a comparison of algorithms developed in [17] and in 
this paper, which we call Alg.l and Alg.2 respectively. 

1. Both algorithms use the nontrivial kernel of the integral operators, but the integral equa- 
tion is inhomogeneous in Alg.l and homogeneous in Alg.2. 

2. Both algorithms use two types of dressing functions: external and internal dressing func- 
tions. However, Alg.l uses two internal dressing functions, $(A, /j) and C(/i;x) with the 
kernel * expressed through $ and C, while Alg.2 uses the single internal dressing function 
\&(A, /i; ), the kernel of the integral operator. 

3. Both algorithms use two disjoint domains on the spectral parameter plane: a continuous 
T> and a discrete D. But the number of points M in D is minimized in Alg.2, while this 
number is an arbitrary M > dimker\l/ in Alg.l. 

4. Alg.2 describes the full solution space of the appropriate n-dimensional nonlinear PDEs, 
while Alg.l describes only an (n — 2) dimensional subspace of analytic solutions. 

5. Both algorithms deal with nonlinear PDEs which may not be considered as necessary 
compatibility condition of some overdetermined linear system of PDEs for the associated 
spectral function. 

6. Both systems of PDEs derived by Alg.l. and Alg.2 have infinitely many commuting flows. 

In this paper we concentrate on the case M — 1. Before giving more details regarding the 
solvability of eqs.(36) and (37), we present the derivation of eq.(l). 
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3 First order quasi-linear PDEs 



3.1 Derivation of the PDEs (1) 

As usual in the dressing philosophy, the x-dependence of the spectral function U (A; x) is in- 
troduced through the x-dependence of the dressing functions. In this paper, as well as in the 
dressing algorithm introduced in [17], we have two types of dressing functions. The internal 
dressing function \I/(A, p;x), appearing in the integral equation (36), and an external dressing 
function G(X; x), A G V U D: 



G(X;x) 



g(X;x), XeV, 
g x (x) A = Zi, 



(45) 



whose role will be explained below. 

This ^-dependence is defined by the equations 



N 



3=1 



N 



G t (X; x) + J2 G Xj {u- x) * A®{v, A) = -TG(v; x) * A^(u; A), 

3=1 

supplemented by the generalized commutation relation: 

A {j) (X,u) * V(is,p;x) = *(A, v\ x) * A (j) (u,p), 



(46) 
(47) 

(48) 



where T is an arbitrary constant matrix. Since A^* and satisfy the same relation as A* 
and *A, it follows that they are expressed as functions of the operators A* and *A respectively: 



A U) * = p ij) (A)*, * A (f} = *p( j \A), 
where p^(-) are scalar functions representable as positive power series: 



(49) 



P°' ) (y) = E c ?V 



fc=0 



*A® = * p V)(A) = Y, 4 * A* • • • *A 

^~ ^ k 

A®* = p (j \A)* = £ c^A*_--*A* 



(50) 



k=0 



From the above definitions and from (44) it follows that 



f a^(X,p), X,peV, 



A U) (X,p) 



A3) 



(A), XeV, p = h, 
«ioV)> A = Zi, peV, 
4i\ A = Zi, p — l\. 



(51) 



Applying the operator A* to eq.(36) and using (37), one shows that A*U is another solution 
of the integral equation (36), i.e.: 



(A* U) = 0. 



(52) 
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Applying instead the operator (d t + J2^ =1 A^d Xi *) to eq.(36) and using the generalized com- 
mutation relation (37) and equation (46), one obtains a third solution of the integral equation 
(36): 

N 

*(\, f ji;x)*LU( f i;x)=0, LU(/i; x) = U t {n; x) + aU) (^ ") * U ^ x )- ( 53 ) 

3=1 

In this section we assume that the three solutions U, A*U and LU of the integral equation 
(36) belong to the same one dimensional matrix subspace spanned by U; i.e.: 

A(X, v) *U {v\ x) = U (A; x)F(x) (54) 

N 

U t (X; x) + J2 ^\X, v) * U Xj [v- x) = U(X; x)F(x), (55) 

3=1 

where the matrices F and F do not depend on the spectral parameters. 

As in the dressing scheme introduced in [17], in order to fix F and F, we use the external 
dressing function G, together with the additional relation: 

G(X; x) * U(X; x) = I, (56) 

where I is the identity matrix. 

Applying G* to equation (54) and using (56), one obtains that 

F(x) = w(x), (57) 

where 

w(x) = G(X; x) * A(X, fj)*U (//; x); (58) 

so that 

A(X,v)*U(v;x) =U(X;x)w(x). (59) 
In addition, applying repeatedly A* to equation (59), we obtain 

p(A) * U = Up(w), (60) 
where p : K. — > R is any scalar analytic function. Therefore equation (55) becomes: 

Ut(X;x) + Ef =1 {U(\;x)pW(w)) Xj = U(X;x)F(x). (61) 
Applying now G* to (61) and using (56) and (60), one obtains 

N 

F(x) = Y,{p ij) W) X] +Tp (0) (wy, (62) 

3=1 
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so that the system system (54), (55) takes the final form: 



A(\,v) * U(v;x) — U(X;x)w(x), 

Ut(\; x) + Ef =1 U Xj (A; x)pW ( w ) = U(\; x)TpP> (w). [M) 

It is a system of overdetermined linear equations for the spectral function U but, as we shall 
see in Sec. 6, its role is different from that played by the usual Lax pair for soliton equations. 
The main feature of the linear equation (63a) is that it does not involve x-derivatives. 
Applying G * A* to (63b), one finally obtains the matrix equation 

N 

Wt + Y, w ^P (j) H = K Tp {0) (w)} (64) 

3=1 

reported in the introduction as equation (1). We recall that w may be either a scalar or a 
matrix. 



3.2 Basic properties of eq.(64) 

The matrix first order quasilinear PDE (64), isolated by the above dressing construction, pos- 
sesses important properties and can be integrated by simple spectral means [21]. Here we 
summarize, for completeness, some of these properties. 

1. We first observe that, if w solves eq.(64), then w T solves the transposed equation 

N 

^ T + P°W)£<. = [p {0 \w T )T T ,w T ]. (65) 
i=i 

where the order of matrix multiplication is inverted. We also observe that T can be chosen to be 
diagonal without loss of generality; indeed, if it were not diagonal, the diagonalizing similarity 
transformation for T would leave eq.(64) invariant, just transforming w according to the same 
similarity transformation. 

2. Let T be the diagonal matrix, and % — 1, . . . , Q be the eigenvalues and the right 
eigenvectors of the matrix w, suitably normalized by the conditions = 1, i = 1,...,Q; 
introduce the matrix E of the eigenvalues: E = diag (e^, . . . , e^) , and the matrix V having 
the eigenvectors as columns (V^- = v\^). Then, if w evolves according to eq.(64) the 
spectrum of w evolves in the following simple way: 

TV 

el fc) + E e SV j) (e (fc) ) = 0, k=l,...,Q, (66) 

N 

V t + J2V Xj p^(E) = [V,T]pM(E). (67) 

3=1 

It follows from (66) that each eigenvalue of w evolves separately according to the scalar 
version of (64), and is constant along the characteristic straight lines 

Xj = pV\eW)t + r h , J = l,..,n, (68) 
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and rjj are arbitrary integration constants. Therefore it is defined by the implicit equation 



(r ?1 ,.., W ) = ( Xl -pW(eW)t,..,x N -pW(e< fc >)f) , 



(69) 



where : JBL n — > R, k — 1, . . . , Q are arbitrary function. Once the eigenvalues are constructed, 
then V satisfies the linear equation (67), with coefficients depending on the corresponding 
eigenvalues. 

3. First reduction. From equation (66) it follows that the condition 

e®( x ) = const, i = l,...,Q, (70) 

represents a symmetry reduction of the system (64). In this case, eq.(67) becomes a linear 
system of (Q 2 — Q) equations with constant coefficients; therefore the transformation from w 
to its spectrum linearizes the flow. 

4. Second reduction. Consider the subspace of QxQ matrices spanned by the basis {u>o, ..,cuq-i} 
given by: 

= I, = ■■ , (uk)ij = 8i+k,j, ■■ , (^Jv-Oy = Si+N-xj, modQ. (71) 

This subspace is left invariant under matrix multiplication, since: 

UjUJk = u k L}j = u j+k , modQ; (72) 

therefore it defines a reduction of (64) from the Q 2 components of w to the Q scalar coefficients 
Vk, k — 1, .., Q of the expansion: 



Q 



w 



(73) 



fc=i 

In particular, if Q — 2, N — 1, p^\y) = y and p^ = 0, equation (64) reduces to the following 
particular case of the gas dynamics equations [20] 



V 2t + V 2 V ixi + V X V 2 



0. 



(74) 



5. The matrix equation (65) can be written as a vector system of the following form 



JV 



w t + ^C«(w)w a;i = B(w) 



(75) 



i=i 



for the Q 2 -dimensional vector w, where 



w 



Wi 



,w 



(Wi---Wq), Wi 



, CW(W) 



B(w) 



Bi(w) 



Bq(w) 



, B fc (w) 



b W (» r )T> T ] fcl 
[^K)T> r ], Q 



Pn 1 



PQI 1 



and 7 is the QxQ identity matrix. By construction, the e^'s are eigenvalues of all matrices CW 
with multiplicity Q. Thus the method of integration of vector first order quasi-linear equations 
developed in [22, 23, 24] may not be applied, at least in the form presented in these references. 
Hereafter T is diagonal constant matrix. 
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4 The general solution of eq.(64) 



In this section we construct, using the dressing method introduced in this paper, the general 
solution of the matrix equation (64), which turns out to be characterized by a nonlinear system 
of non-differential equations for the components of matrix w in the following way. 

Proposition 1. Let : M. N — > R, i, j = 1, . . . , Q be Q 2 arbitrary scalar functions, repre- 
sentable by positive power series, so that F i j(M 1 , . . . ,M N ) are well defined matrix functions, 
where Mi, ... , M N are arbitrary Q x Q matrices. Let . . . , Tq} be the elements of the con- 
stant diagonal matrix T. Then the general solution of the matrix PDE (64) is characterized by 
the following system of 2Q 2 non-differential equations: 

Q Q 

w a( 3 = Y] Y] {(ui^y^aslF&nixiI - p (1 \w)t,...,x N I - p (N) (w)t)) (ui(z)) 7172 , (76) 

. . \ ' 72/3 



V(«i(^)U(e^ (0) ^) =5 aPl a,P = l,...,Q. (77) 

for the components of the matrix solution w(x) and of the auxiliary matrix function u\(x). 
Remarks 

1. If T = 0, equation (77) gives u\ = I; then equation (76) reduces to the following system of 
Q 2 non-differential equations for the components of w: 

Q 

w a p = J2(F ai (x 1 I-p( 1 \w)t,...,x N I-p( N \w)t)) , a,P = l,...,Q. (78) 

7=1 

This equation can also be constructed using equations (66) and (67) for T = [21]. 

3. If, at last, one is interested in the scalar version of equation (64), its general solution, 
characterized by the scalar version of (78): 

w = f(xi — p {1) (w)t, ...,x N - p (N) (w)tj , (79) 

is also easily obtainable using the method of characteristics. 

4. Due to the presence of the Q 2 arbitrary scalar functions Fij, i, j — 1 . . . , Q, the above non- 
differential equations characterize the general solution of the matrix PDE (64). In particular, 
if one is interested in solving the Cauchy problem in ~R N with the prescribed initial condition 
w (xi, . . . = w(xi, . . . ,xn,0), equation (77), evaluated at t — 0, implies that Ui\ t=0 = I. 
Then equation (76) at t — implies that 

F(x 1 ,...,x N )=w (x 1 ,...,x N ). (80) 

Known the functions F^, the nonlinear non-differential equations (76,77) allow one to construct, 
V£, the solution w(xi, . . . , Xn, t). 
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4.1 Proof of Proposition 1 

To prove Proposition 1, we make use of the main ingredients of the dressing scheme: eqs.(54),(57) 
and (56), written explicitly using (41,44,45): 

a(\, fi) * u(fi; x) + aoi(\)ui(x) = u(\; x)w(x), A G V, (81) 
a io(/-0 * u (fi] x ) + dnUi(x) = Ui(x)w(x), A = li, 

g(li;x)*u(Li;x)+g 1 (x)u 1 (x) = I, (82) 

eq.(58) for w: 

w(x) = G(X;x) * A(X, /i) * U(/i;x) — g(X;x) * a(X, /i) * u(/i;x) + (83) 
g(X; x) * a i(A)ui(a;) + gi{x) (a w (X) * u(X; x) + an«i(a;)j , 

and the generalized commutation relation (37), which we rewrite to emphasize the fact that 
the parameters A and \i take values in different domains: 

A{X, v) * ii) — *(A, v) * A{v, //), XeV, lieVUD. (84) 

Using formulae (40), (41), (43) and (44), eq.(84) takes the form: 

a(A, v) * ip(is, il\ x) — ip(X, v\ x) * a(v, li) + ipoi(X; x)a w (Li), X,fxEV, (85) 
a(A, v) * t/; 01 (u; x) = ip(X, v\ x) * a m (v) + ^oi(A; x)a n , XeV, /i = l ± . (86) 

This system should be viewed as a system of linear equations for the functions ip and > where 
largely arbitrary. The following two choices for a and a have been explored so far: 

1. a(A, li) — a(X)8(X — /i), a(X, li) — a(X)5(X — [/), (87) 

2. a(A,/i) = iS'(X -/i), a(X,Li) = -i6'(X- li), (88) 

where 5 and 8' are the Dirac function and its derivative. 

It is possible to verify that the first choice (87) corresponds to the case in which the dressing 
function ^ does not depend on x and, consequently, the eigenvalues of w are constant as well. 
Therefore the choice (87) leads to the trivial reduction discussed in the previous section (see 
(70). 

As we shall see in the following, the second choice (88) allows one to capture instead the 
general solution of the matrix PDE (64). 

4.1.1 Solution space associated with eqs.(88) 

Construction of the dressing function \I/(A, fi;x) Under the assumption (88), function \& 
is completely defined by the equation (46), which reads 

N 

iPt(X,Li;x) + J2p U) ( id x)^ j (X,^,x)^0, X,liEV, (89) 
j'=i 

N 

ii ou (X;x) + J2p U) (iWoi Xj (X;x) = 0, X E V, lc = h, (90) 
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and by the system (85,86): 

itpx(X,i^;x) = itp^/i-x) + tp 01 (\;x)a w (ij), X,fieV, (91) 
ii[> 01x (X;x) =ijj(\,v;x) * a 01 (v) + ip 01 (X;x)a u , XeV, n = h. (92) 

Equations (89-92) suggest to represent ip, ipoi and a w in the Fourier forms 

OO OO N 

dx l dq I dk tfj(x,q,k)e j=1 , (93) 



-OO — OO R« 

OO N 



^bi(A;x)= dx dkip 01 (x,k)e j=1 , (94) 



-OO 



aM = / a w (q)e iq »dq, (95) 



-OO 
OO 



u(X;x)= J ^u(q;x)e iq \ (96) 

— OO 

where k = (ki, . . . , k^). 

For future convenience, we also represent a i(/i) via the contour integral 

aoM = 7?~ I aoi(q)e iq »dq, (97) 



2ni 

r i 



where r i is a sufficiently large contour containing all singularities of the integrand. 

With the above representations, eqs. (89,90) are automatically satisfied, while eqs. (91,92) 
yield the relations: 

tjj(x, q, k)(x - q) = -Tp 01 (x, k)d 10 (q), (98) 

OO OO 

tp 01 (x,k)(x + a u ) = - J dfi J dqij(x,q,k)e iq < 1 a 01 ( y {i). (99) 



— oo — OO 



We solve the first of these equations with respect to ip 

ijj(x,q,k) = + (f)(x,k)5(x-q), (100) 

x — q 

where is an arbitrary function of its arguments, and substitute the result in (99): 

OO 

4> 01 (x,k)(x+a u -?7(x)) = - J ^0(x,/c)e^a O i(/i), (101) 



— OO 
OO OO 



rj(x) = J dji J dq 



a w (q)e iq ^a 01 (iJ,) 
x — q 



-oo — OO 



16 



We take the following solution of this equation: 



— oo 

fi(x) = x + an — 7)(x); 



(102) 



then eq.(100) yields 

Tp(x,q, k) = 0(x, k) 



e^ a „ lMn -'(x)a lo(g ) + 
x — g 



(103) 



Substitute ^ and V'oi into (42), with M = 1, and apply the operator J d\ e'^ x - to the result, 

— oo 

obtaining: 



*(-t,x 1 -pW(t)t,...,x N -pW(£) t } x (104) 



u 



where 



0, 



^(-e^i-p (1) (0*,---,^-p (jV) (0*) = / ^-^^e^^^^-^^dfc. (105) 

R N 

We also remark that eq.(81) takes the form 

— iu\(\; x) + aoi(\)ui(x) — u(X; x)w(x), A e T>, (106) 
oio(aO * x ) + ttnUi(x) = A = l\. (107) 

Due to equation (97), the solution of (106) leads to the following representation of u: 

u(X;x) = ^-j dqa 01 {q)e iqX Ul (x)(w(x) - qij \ (108) 
r i 

different from the Fourier representation (96). Then eq.(107) yields: 

oo 

- — ; J dfi <j> dq a w (n)d i(q)e iqfJ 'Ui(x) (^w(x) — q j = u\(x)w(x) — a u Ui{x). (109) 



-oo Toi 



The formulae derived so far are applicable to both the scalar and the matrix equations. Below 
we consider these two cases separately, starting with the simpler one. 
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Scalar nonlinear PDEs In the scalar version 

N 

w t + ^w Xj pW(w) = (110) 
j'=i 

of (64) the RHS is zero. In this case we can set T = 0; therefore equation (47) admits the 
trivial solution g — 0, g± — 1; so that u± — 1 (see (82)). 

It follows that eq.(83) coincides with eq.(107), while equation (109) reads: 

oo 

- \dAdg aM ~ aMeim = W - a n . (Ill) 
2m J J w-q 



-oo r i 



This equation indicates that, if a w and a i are given, then (111) is an algebraic constraint for 
w, which would imply that w is constant. To avoid such a trivialization, eq.(lll) must be taken 
as a definition of aio or aoi, and w must be considered as the independent variable. 

In addition, one should make sure that the condition (111) does not coincide with the 
condition Q(x) = (see (102b)); i.e.: 

oo oo 

/ x_ f . f > a 10 (q)e iqfM a 01 ([i) 

r](x) = dfi dq f x+flu, (112) 

J J K-q 

— oo — oo 

otherwise equations (102)- (104) would make no sense. 
A possible (and simple) choice for a m is given by: 

«oi(9) = — ~~T> (H3) 
q — o 

implying 

a M = e i6 ". (114) 
Then, assuming that aio(o) is an entire function, eq.(lll) gives 

- ( \ (w + a u )(w + b) 

aio{w) = haio(-o), (115) 

Z7T 



while eq.(112) reads 



??(x) = 27r Q ^Y^ + Q(x)^0. (116) 



Hereafter we assume, without loss of generality, that 

a l0 (-b) = =^ r](x) = 0, fi(x) = x + a n . (117) 
Now we rewrite eq.(104) using eqs. (111-115): 

( - £, *! - P « (0*, - P w [«(£; *) 

oo \ 
-oo / 



118) 
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Using again equations (114), eq.(108) yields u in terms of w: 



so that 



p ib\ iwX 

u(X;x) = — , (119) 

w — b 



u{ q]X ) = 2 ^ q - h) - 5 [ q - W \ (120) 
w — b 

Substituting it into the eq.(118) one gets: 

( - e, *i - p (1) (o*, ...,x N - P ™ (ot) = o, (i2i) 

which is satisfied iff 

^-tu.^"/'^....!^^^)*) =0. (122) 

The implicit eq.(122) for w suggests to take function (p(y,Xi, . . . , Xjv) i n the form: 

(p(y, x u ...,x N ) = y + F(xi, ...,x N ). (123) 

Thus eq.(122) yields 

w = F(x 1 - p (1) (w)t, ...,x N - p {N) (w)tj , (124) 

which is the well-known non-differential equation defining implicitly the solution of the Cauchy 
problem in 

N 

w t + y^ j w x .p^ J \w) = 0, w\ t=0 = F(x 1 ,...,x N ) (125) 

3=1 

for the scalar version of equation (64). 

Then the direct problem, the mapping from w\ t=0 to the dressing function ip (or 0, via (102) 
and (103)), is simply given by 

w\ t=0 = F(xi, . . . ,x N ) =>• <f>(y,x 1 ,...,x N )=y + F(xi,...,x N ). (126) 

From the inverse problem point of view, given from the initial data through eqs. (122,123), 
the spectral function u is reconstructed solving equation (118), which is equivalent to 

oo 

*K; *) - I ^ ''- t) ''- t "' % *) = - t) + - «,i27) 

7 2ir [a u - £J(g - an - £ 



— oo 



where a(x) is found requiring that u be compatible with the expression (120). This request 
fixes the value a = —2n/ (w — b). 
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Matrix nonlinear PDEs In the matrix case, we choose the operators A* and *A to be 
scalar; i.e., a 10 (X) and a i(A) are scalar functions and a u is a scalar parameter. Now we 
consider the general case Ui ^ I. Then eq.(108) yields 

u(\; x) = j dq a 01 (q) (w(x) - q^ ' e iqX . (128) 

r i 

Assuming the invertibility of u 1: eq.(109) yields the following matrix generalization of (111): 



2m 



oo 

J dfi j) dq ai (fi)a i(q) e im (w (x) — g/) =w(x) — a 11 I. (129) 



-oo r ( )i 



The expressions for t], doi(q) and aoi(p) remain the same, see eqs. (112-117). 

We also assume that w is diagonalizable with eigenvalues {e^\ . . . , e^}, and that V(x) is 
the matrix of right eigenvectors, so that: 

w = \/diag(e (1) ,...,e (0) )\/- 1 . (130) 

Then, applying V~ l and V respectively from the left and from the right to equation (129), we 
obtain the scalar equation (111), in which w is replaced by its eigenvalues. Consequently, also 
the expression of 5io is given by (115), and (117) holds too. 

The matrix generalizations of equations (118), (119) and (120) read, respectively, 

P (1) (£)t,...,x N - P W (Ot) fe; x) - 

00 X 1 (131) 



u(fi; x) = Ul (x)V(x) (e^ 7 - e^ s ) (E - bl^V'^x) = Ul (x)(e ib ^ - e^ w )(w - bl)' 1 , (132) 
and 

u(q; x) = 2n Ul (x)V(x) (j5(q - b) - 5{ql - £)) (E - bl^V'^x) = 
2nu 1 (x) (j5(q -b)- 5(ql - iu)) (w - bl)~ l . 

Substituting (133) into the integral equation (131) and multiplying the result by the matrix 
(w — bl) from the right, one obtains the matrix distribution equation: 

4> ( - Z,X! - p« , • • • , xn ~ P iN) (0*) u 1 (x)5(£I-w) = 0. (134) 
Since 5(£I — w) = V5(£I — E^' 1 , equation (134) is equivalent to the distribution equation 

£? 1= i £2=i E? 3 =i *«n ( - 6 *i - fF> (Ot, ...,x N - (Ot) x 

V J (135) 

( Ml (x)) 7l72 \/ 7273 (x)5(e - e^)V~l{x) = 0, 
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which implies 



(136) 



^71=1 ^72=1 ^73=1 ^^71 

(Mi(x)) 7l72 \/ 7273 (a;)K 7 ^(x) =0, a,P=l,...,Q. 

Choosing, in analogy with the scalar case, the following form for the matrix <fr: 

fe(j/,a;i,...,a;jv) = yhj + F i:j (x u . . . , x N ^j , (137) 

and eliminating eigenvalues and eigenvectors, equation (136) becomes 
Q Q 

{uiw) aP = V Y^Mx^jF^ixJ - p (1 \w)t, . . . ,x N I - p (N \w)t)) . (138) 



71=1 72=1 



Applying w : 1 from the left we finally obtain the non-differential equation (76), involving the 
solution w(x) of (64) and the auxiliary matrix function ui(x). 

4.1.2 The dressing function G(X;x) and the construction of u±(x) 

If T ^ 0, the solution of (64) is computed involving also functions U\ and G, which are related 
by equation (82). Substituting u given by (132) into (82) and using the following Fourier 
representation for g{ji; x): 

oo 

/J 
■£g(x;x)e 1 ^, (139) 



one obtains 



Yl Yl [{9 ail (-b;x)I - g ail (-w]x)^j(w - W) -1 ] (uifa))^ + (140) 

7i =1 72=l 

(gi(x)ui(x)) a p = Sap, a, (3 = 1, . . . ,Q. 

This is a determined system of linear equations for the elements of U\, where functions g and 
g 1 are solutions of eq.(47), which reads, in terms of (44,45) and (51): 

N N 

g t (\;x) +Y.9*Mi x ) * aU) ("> A ) +I>i*>)4 ) (A) = (141) 

3=1 3=1 

-Tg(v- x) * a^{v, X) - T 9l {x)a^{\), XeV, 

N N 

9it(x) +Y,9-M X ) * 4V) = (142) 

3=1 3=1 

-Tg^af} - Tg(fj,; x) * a£V)» X = h- 

The representations of functions j = 0, 1,2, .. in terms of A can be obtained recursively 
through the formulae 

a (ri) (A, v) * a(u,fi) + a$ (X)a w (fi), X, fi eV, 



A (n) (X,v)*A(v, n) 



a { ^ ) {X)a 11 + a^(X,iy) *a 01 (v), XeV, ji = I 



a^(u) * a(u,fx) + a$a m ((x), X = h, fxeV, 

k a$(v) * a 01 (u) + a$a u , X = fi = l v 



(143) 
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Due to the choice (88) for a, and to the representations (95,113,114) of a±o and aoi, it is possible 
to look for a solution of the system (141, 142) satisfying the following properties: 



0i = O, g(\;x)*a$(\) = 0, 

g(n; x) * a( j \n, A) = g(p; x) * p U) {a{n, A)) = pW(id x )g(\; x), j = 1,2,... 
which reduce such system to the single equation 



(144) 



N 

g t (X;x) +J2p U) (idx) g Xj (\;x) +Tp^(id x ) g(X;x) = 0. (145) 
j'=i 

Indeed, as it is possible to verify using equation a 01 (p) = exp(ibp) and formulae (143), the 
restrictions (144) are simply satisfied by imposing the single condition g(—b;x) = 0. 

Then the matrix function g satisfying (144,145) can be represented in the following Fourier 
form: 

oo 

g(X- x ) = J g J dk J duo g{x,k,u)e Mi ^=^- iujt , 

-oo w rn r w (146) 

g aP (x, k, oo) = (x + b)h a p(k)5(uo - Ef =1 k^\-H) + ip(°)(-x)T Q ) 

where T u is any contour passing through the support of the Dirac function, and h a @, a, (5 — 
1, . . . , Q are arbitrary scalar functions of their arguments. 
Then replacing 

g{x;x) = J dk J doo ~g{K,k,uo)e^^ k ^- iult . (147) 

in equation (140), one obtains: 

Q 

V \H ail (x 1 I-p^\w)t,...,x N I-^ N \w)t)e-^ ) ^M K(x)) 7172 = 5 aP , (148) 



71,72=1 

where 



H afS ( Xl ,...,x N )= [ dkh a/3 (k)e l ^ klXl . (149) 

If we choose h a p{^, k) = (f| ■ 5(kj))5 a p, =>• -f^ Q /3 = <W; then equation (148) reduces to (77). 
This simplification can be done without loss of generality since, if it were not made, it would 
correspond to the following redefinition of the arbitrary matrix function F appearing in (76): 

F -> HFH~\ □ (150) 
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5 Second order quasilinear PDEs 

5.1 Derivation of second order PDEs 

To increase the order of the nonlinear PDEs, we introduce the x-dependence through the next 
equations 

N 

%J\V,x) + J2^ mj \X,iy) * **>,A*;aO = °, ( 151 ) 
3=1 

N 

G t jX, q;x)+Y, G Xj (v, q; x) * A™(v, A) = -q m S^G(u, q; x) * A^°\u, A), (152) 
3=1 

supplemented by the generalized commutation relation: 

A {mj \\v) *V{y,n;x) = V(X,is;x)*A( mj \is,p J ), (153) 

where q is a new spectral parameter, are constant diagonal matrices, 

A (mj h = p {mj \A)*, *A {mj) = *p {mj \A), (154) 

and p( m ^ : R — > R are arbitrary analytic functions. 

Together with the field w introduced in the previous section, we introduce also the fields 

y( n \x) — (G(/jl, q; x)q n ) * U(fi; x), n = l,2,..., v(x) = v (1) (x), (155) 

where now "*" means integration also over q. 

Applying the operators (d tm + Yl!j=i A ( - m ^d Xj * j to (36), one gets 

N 

M>(A, p) * (U tm (p; x) + Y, ^ {m3) (», v) * U Xj (y, p; x)) = 0. (156) 

3=1 

Assuming, as before, that the solutions of (36) belong to the one-dimensional matrix subspace 
generated by U, the spectral equations are similar to (54,55): 

A(X, v) * U(u; x) — U(X; x)F(x), (157) 

N 

Ut m (A; x) + A(mj) (A, v) * U Xj (u, p; x) = U(X; x)F™ (x) . (158) 

3=1 

Applying G* to (157) and (158) and using the condition (56) and equations (152), (60), we 
obtain the expression of F(x) and F^ m \x) in terms of the matrix fields w and v^: 

N 

F(x) = w(x), F {m \x) = {P {mj \w)) x . + S {m) v {m) {x)p {m0) {w). (159) 

3=1 

Therefore equations (157) and (158) become the following overdetermined system for the spec- 
tral function U(X;x): 

A(X, p)U(p;x) = U(X;x)w(x), 

Ut m {X ]X ) + J2 U Xj {X-x)p^\w) = U(X;x)S (m ^v^(x)p^(w). (160) 

3=1 
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In order to construct the nonlinear PDE for w, we apply G * A* to equations (160b) and use 
equations (60), (151) and (155), obtaining 

N 

wt m + J2 w ^p (m3) ( w ) = Ks (m V"V m °)H]. (i6i) 

3=1 

To write the equations for v^ k \ we apply instead (Gq k )* to eq.(160b), obtaining: 

N 

vi k J +J2 v x-P imj) ( w ) + S {m) v ik+m) p im0 \w) = v ik) S {m) v {m) p im0) {w). (162) 

3=1 

To construct a complete system of nonlinear PDEs, consider the three equations: 

N 



(w) = V WsW v WpM(w), (163) 

3=1 
N 

3=1 
N 

+ y^fffp (2j) (w) + s ,(2) w (3) p (20) ( w ) = v 

(w), 

3=1 

corresponding to (m, k) = (1,1), (1, 2), (2, 1). 

We may choose, without loss of generality, = 1, = S. Then, expressing and 
v ^ in terms of v — v^: 

vW=v*-CM{P {W) {™))-\ nfi4 x 

V (3) = V (% - A^Xp^H)- 1 , ^ 

one obtains a single matrix equation of the second order for v and w: 



C 2 (v)p (10 \w) + SA(A(^)(p (10 H™)) _1 JP (20 H™) = SC^p^iw) + [SC^p^iw)^] + 
Sd(v) [(p^(w))- 1 v,p^(w)p^(w)] + [v,Sv 2 ]pW(w)p (10 \w), 

(165) 

where the differential operators C m , m = 1, 2 are defined by 



N 



C m (f(x)) = f t Jx) +J2Ux)P {mi) W- (166) 



i=l 



A complete system of second order matrix PDEs in arbitrary dimensions for the fields w 
and v is obtained coupling equation (165) with equation (161) for m — 2, obtaining: 

C 2 (w) = [w, S(v 2 - C^Mp^iw))- 1 ], 

C 2 (v)pW(w) + SCi (Ci (v) (pW (w))' 1 ^ pW (w) = SXi^V^H + [SC^p^iw)^] + 
Sd(v) [(p^(w))~ 1 v,p^(w)p^(w)] + [v,Sv 2 ]pW(w)p (10 \w). 

(167) 

But the independent equation (161) for m = 1: 

d(w) = {w,vp {10) {w)} (168) 
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is also satisfied by the fields w and v and must be viewed as an integrable constraint for the 
evolutionary (with respect to the time t 2 ) system (167). 

We remark that the evolutionary system in arbitrary dimensions (167) is completely inte- 
grable only under the nonlinear constraint (168). 

Here we consider two explicit reductions. 

1. Let 

N = 2, p M=pW=w, pW=p^=I, pW = p(2D = o, (169) 

then 

&(f) = ft, + f Xl w, C 2 (f) = f t2 + f X2 w, (170) 

and one obtains the nonlinear system: 

w t2 + w X2 w = [w,S(v 2 - v tl - v Xl w)], (171) 
v t2 + v X2 w + S(v tl + v Xl w) tl + S(v tl + v Xl w) Xl w = 

s {iv 2 )ti + (v 2 ) Xl w^j + [S(v tl +v xi w),v] + [v, Sv 2 ], 
subjected to the constraint 

w h +w Xl w = [w,v\. (172) 
Using this constraint, eq.(171b) can be rewritten in the more convenient form 

v t2 + Sv tltl + v X2 w + Sv xixi w 2 + 2Sv xltl w - [S,v](v tl + v xi w) - (173) 

2S(v tl v + v Xl vw) + [Sv 2 , v] = 0. 

2. Let 

N = 2, p^ = p^ = p< 10 > = p(20) = 7> p (12) =p (21) =0> (174) 

then 

Ci(f) = f tl + f Xl = /n, = ft 2 + f X2 = f T2 , (175) 

and the nonlinear system reads 

w T2 = [w, S(v 2 - v n )], (176) 
v T2 + Sv T1T1 = S(v 2 ) T1 + [Sv T1 ,v] + [v, Sv 2 }, 

with the constraint 

w Tl = [w,v]. (177) 

We remark that, in this reduction, equation (176b), which involves only the field v, is the matrix 
Burgers equation [26]. 

Using equations (41), (45) and (155), the fields have the following spectral representa- 
tion: 

^) = (G^)*V M = + / 9l(5 ,,) 9 M 9M , (178) 



25 



The dressing functions ip and G are responsible for the dimensionality of the space of 
solutions for w and v respectively. 

As we did in Sec. (3.2), we may separate diagonal and off-diagonal parts of the equation 
(161), obtaining 

N 

e£ ) + E e SV ww " ) (e (fc) ) = 0, k = l,...,Q, (179) 
j'=i 

N 

V tm + V X] p [m3 \E) + S (m) v {rn) V = VD, D diagonal. (180) 

3=1 

Equation (179) coincides with equation (66); i.e., the eigenvalues of w evolve separately accord- 
ing to the scalar version of equation (161). Instead, eq. (180) cannot be written in the same 
form of eq. (67), since v ^ is not diagonal matrix. This is a principal difference between eq. 
(161) and eq. (64). 



5.2 The solution of the system (167), (168) 

In this section we construct, using the dressing method introduced in this paper, the solution 
of the matrix equations (168) and (167), which turns out to be characterized by a nonlinear 
system of non-differential equations for the components of the unknown matrices w and v in 
the following way. 

Proposition 2. Let F {j : R N -> R, i, j = 1, . . . , Q and Hij : R N+1 — > R, %, j = 1, . . . , Q be 
2Q 2 arbitrary scalar functions, representable by positive power series, so that Fy(Mi, . . . , M N ) 
and Hij(q, Mi, . . . , Mn) are well defined matrix functions, where Mi, . . . , are arbitrary 
Q x Q matrices. Then the solutions of the matrix system (167), subjected to the matrix 
constraint (168), are characterized by the following system of 3Q 2 non-differential equations: 

«v = Y? n ^ 5= ^i\x)) aS (F Sll (xii - ELi p (ml) H*m, ...,x N i- ELiP (mA °MO) 72/3 x 

(ui(x)) 1112 = 0, a, (3 = l,...,Q, 

(181) 

OO p 

/ dq E? 1= iE? 2= i ^ 71 (g^i/-ELiP (ml) H^,---,^-ELiP (mJV) HUx 

— oo 

e-^ s "'^ Mw H ( Ul (x)) 7lT2 = 8*0, a,(3=l,...,Q, 

-I 72/3 

(182) 

oo . 

Vafi(x) = j dqqJ2l^ =1 {H ail (q, xj - ELi /> M (*- • • • ^nI ~ ELi P {mN) W^) x 

— oo 

/ 72/3 

(183) 

for the components of the unknown matrices w(x) and v(x), and of the auxiliary matrix function 
ui(x). 

Proof of Proposition 2. To prove this Proposition, we proceed as in Sees. 3,4. The solutions 
w and v^> are constructed using the algorithm presented in Sec. 4, but the functions ^> and G 
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will be now denned by eqs. (151,152). Eq.(151) yields 

N 

^ t J\,l,;x) + J2p {mj) (i9 x )^ Xj (X,i,;x)^0, XeV, (184) 

3=1 
N 

Voi tm (A;x) +J2p {mj) (id x )^ 01x .(X;x) = 0, A = h, 

3=1 

where m = 1,2. The analysis of the system (184) coincides with the analysis carried out in 
Sec. 4. 1.1. but ip and ^ i exhibit the following Fourier representations: 

oo oo / \ 

ip{X,H;x)= / dx dq dkip{x,q,k)e \ /, (185) 



— oo — oo 
oo 



i/VA;x) = dx dk ifj 01 (x,k)e V /, (186) 

-OO RJV 

and equation (76) is replaced by (181). 

The x-dependence of G is introduced by eq.(152), with a(A, p) given in (88): 

AT N 

g tm (A, g; x) + £ ^ (i/, g; x) * (i/, A) + ]T <^ (g, x)a^ } ( A) = (187) 

3=1 3=1 
_ S (m) q m ^ g . x) _ a(fn0)(l/j A) + ^ X ) a (-°)( A )) , \eV 

N N 
3=1 3=1 

-S^q m (gi(q,x)a^ 0) +g(»,q;x) *4™ 0) (/i)), A = h, 
where m = 1,2. We impose again the conditions 

gi (q;x)=0, g(fi.,q;x) *a£ lj) (/x) = 0, (188) 
g(p, g; x) * a^)(//, A) = </(//, g; x) * (a(/i, A)) = ^'>(<9 A )»(A, g; x), j = 1, 2, . . . 

on the solution of the system (187), reducing it to the single equation 

N 

9tJ\ q- x) +J2p {mj) (idx) g xj (\, q; x) + q m p imQ \id x ) g(X, q; x) = 0. (189) 

3=1 

Functions g satisfying (188,189) can be written in the following Fourier form: 



OO oo 



g(\,q;x) = J dk J du J dx g(x, q, k, uj)e i>cX+i '^'= 1 t|Ir '^i u(mlim (190) 



]RJV — oo — oo 

where 

N 



g a p(x,q,k,u;) = (x + b)h aP (q,k) J] sfa™ - £ k lP ^ l \-x) + iS^q m p^\-x)). (191) 



m=l 1=1 
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Now g(x,q] x) takes the form: 

oo 

g(x,q;x) = J dk J du g(x, q, k, u)e ij: ^ ^r^Li^- 1 '^ (192) 
i" -00 

and the analogue of equation (77) is equation (182), where 

H{q,x u ...,x N ) = J dk h(q, k)e ij: ^ k ^ . (193) 

At last, substituting in equation (178) the expressions (132) and (190) of u and g, we obtain 
(183). □ 

Initial-Boundary Value Problem. A well-posed initial-boundary value problem for equa- 
tions (167), (168) is the Cauchy problem for the system (167), in which the initial conditions for 
w and v at t 2 — are any pair of matrix functions satisfying the multidimensional PDE (168). 
This initial constraint can be satisfied assigning arbitrarily v at t 2 = and w at t± — t 2 — 0: 

w m = w ( x )\ ti=t2=Q „(o) = V ( x )\ t2=0 . (194) 

Then, integrating equation (168) with respect to t±, one obtains = w(x)\ t2 =o, V£i. At last, 
given v(°\w(°\ the evolutionary system (167) allows one to construct w and v, \/t 2 . 

The algorithm allowing one to integrate such an initial-boundary value problem consists of 
three steps. 

1. At t 1= t 2 = 0, the system (181,182,183) reads: 

w m = ( u myi F ( Xu . . . , XN )u(°°\ (195) 
H(0,x 1 ,...,x N )u m = I, (196) 
= H'{0,x l ,...,x N )u^ \ (197) 

where = v^(x)\ tl =o is given, = ui(x)\ tl =t 2 =o, and 



H(ti,X!,...,x N )= f dq H(q,X!,. . .,x N )e qtl , 

—00 

H'(0,xi, ...,x N ) = H tl (t 1 ,x 1 , . . .,x N )\ tl=0 . 



(198) 



This system of three matrix equations must be solved for F(xi, . . . , xjy), and H'(0, x±, . . . , xn). 
Function H(0, x 1 , . . . , x N ) remains arbitrary. 

2. The system (181,182,183), evaluated at t 2 = 0, reads: 



Q 

,(o) 

71,72,5=1 



«■:.■';= E (( w(0) ) _1 U ^7 1 (^i^-p (11) (^ (0) )ii ) ---,^-p (1JV) (^ (0) )ii) (« (0) ) 7172 ,(i99) 

= 1, . . .,Q, 

E E [^ T1 (p (10) (^ 0) )ti,xi/-P (11) (^ 0) )ti,..., (200) 
71=1 72=1 

x N I - p^iw^h)} (u (0) ) 7172 = <f a/J , a,/3 = l,...,Q, 
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Q Q 



v $ = EE [^ 71 (P (10) (^ (0) )ti,-i/-P (11) (^ (0) )ti, 



(201) 



71=1 72=1 



x N I - p^iw^t,) (« (0) ) 7172 , a,/3 = l,...,Q, 

-I 72/3 



where = Mi(x)|t 2 =o- Eqs. (200,201) must be solved for and H ail (ti, zi, . . . ,zn). Then 
eq.(199) gives w^°\ 

3. After that, the functions w, U\ and v can be constructed as the solutions of the non- 
differential system (181,182,183), Vt 1? t 2 . 



6 Auxiliary linear system 



We have seen in the previous sections that the dressing algorithm produces, together with the 
integrable nonlinear PDEs in arbitrary dimensions, also the associated linear overdetermined 
system of equations for the spectral function U, whose coefficients are related to the fields of 
the nonlinear PDEs. 

In this section we show that, unlike the classical S-integrable case, the integrability 
condition for the overdetermined system for the spectral function U does not fix 
completely the integrable nonlinear PDEs. This is not surprising, since the derivation 
of the nonlinear PDEs of this paper requires, together with the linear system for U, also the 
external dressing function G and the constraint (56) for it. 

Consider the following linear system, which corresponds to the general form of the nonlinear 
PDEs treated in this paper: 

£ (0) := A(\,n) * U{n\x) — U(X;x)w(x), (202) 

N 

:= U tm {\-x) + Y,U X] {\;x)p^\x) = U{\-x)F^\x) (203) 

3=1 

for the arbitrary matrix coefficients w, and F^ m \ The compatibility of equation (203) 
with (202), for any fixed m, yields: 

N N 

U Xj [w, pW\ + U [w tm + ™*/ m3) - K F {m) }] = 0. (204) 

3=1 3=1 

Assuming that U and its derivatives are independent matrix functions, equation (204) implies 
the following relations among the coefficients of the system (203,202): 

[w,p( mj )] = 0, (205) 

N 

™ tm + w */ mj) - \™> F(m) ] = °- ( 206 ) 

3=1 

In the scalar case, eq.(205) is identically satisfied, so that eq.(206) is some relation among 
p( mj \ w and for any m. Thus, the system (205,206) cannot be considered as a complete 
system of nonlinear equations for some fields. 
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In the matrix case, eq.(205) implies that the part's are arbitrary functions of w, repre- 
sentable in a power series of w with scalar coefficients depending arbitrarily on x. The par- 
ticular case in which these coefficients are independent of x, is in agreement with the form 
of the nonlinear PDEs derived in this paper, but no condition is imposed on the coefficients 
p{m) ^ Therefore, also in this case, eqs. (205,206) cannot be considered as a complete system of 
nonlinear equations. 

The compatibility of two equations from the list (203), for instance, and E^ n \ n ^ m, 
yields instead: 



N N 
i=l j=l 



(207) 



N 



3=1 



N 



E (p. 



(nj)p(mi) _ p{mj) p{ni) 



+ 



U 



N 



F t ] - F t ] + J2( F t ] p {ni) - F ifp {mi) ) + \ p{n \ F 



(m)l 



i=l 



Assuming again the independence of U and its derivatives, one obtains the following relations 
among the coefficients: 



[ p (nj\ p (mi)j = 0, 

P^ + E P^ ] P (mi) ~ [P (nj) , F (m) } = m, 

i=l 

F t + E F x7P (mi) - FWFW = m. 



(208) 



Equation (208a) follows from the relations (205), and is satisfied if, for instance, p*™*-* are 
functions of w, as it has been previously established. Then equation (208b) is a consequence 
of the relation (206). Eq.(208c) prescribes relations among the F^'s, but leaving one of 
them free. Therefore the system (205-206,208) cannot be considered as a complete system of 
nonlinear PDEs for the coefficients of the linear system (202,203). 



7 Conclusions 

Using a new version of the dressing method, based on a homogeneous integral equation with 
nontrivial kernel, we constructed a new type of integrable multidimensional nonlinear PDEs. 
There are no formal restrictions on the dimensionality of the PDEs, while these restrictions 
are very severe in the case of the classical S'-integrable systems, which are known to be the 
first examples of nonlinear PDEs treatable by dressing technics. Several modifications and 
extensions of the dressing algorithm presented here will be considered in subsequent papers. 
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